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The SciBooNE Collaboration reports a measurement of neutral current coherent tt" production 
on carbon by a muon neutrino beam with average energy 0.8 GeV. The separation of coherent from 
inclusive tt" production has been improved by detecting recoil protons from resonant tt" production. 
We measure the ratio of the neutral current coherent 7r° production to total charged current cross 
sections to be (1.16 ± 0.24) x lO^'^. The ratio of charged current coherent n^ to neutral current 
coherent 7r° production is calculated to be 0.14+Q2g, using our published charged current coherent 
pion measurement. 
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I. INTRODUCTION 

Recent measurements of coherent pion production by 
muon neutrinos at neutrino energies around 1 GeV have 
inspired significant discussion |l|. In coherent pion pro- 
duction, the neutrino interacts with an entire nucleus; no 
nucleon recoil occurs and the tt" tends to be emitted in 
the forward direction. 

For charged current (CC) coherent pion production, 



-f-fPrcsent address: Department of Physics, University of Tokyo, 
Tokyo 113-0033, Japan 



both K2K and SciBooNE set limits on the ratio of CC 
coherent pion production to the total CC cross sections 
near 1 GeV[2, Q, These published upper limits are sig- 
nificantly lower than those predicted by the Rein and 
Sehgal model [61, |7| which is widely used for many neu- 
trino oscillation experiments. Meanwhile, evidence for 
neutral current (NC) coherent pion production with neu- 
trino energy less than 2 GcV has been reported by the 
MiniBooNE Collaboration Q. The SciBooNE collabora- 
tion also reported nonzero NC coherent pion production 
Q although the result is only 1.6 standard deviations 
above zero coherent production. Currently there is no 
theoretical model which can accommodate all of these 
recent measurements. Further experimental inputs may 
help the development of theoretical models. 

NC coherent pion production at neutrino energies 
around 1 GeV is also important for neutrino oscillation 
experiments as a substantial contribution to NC tt" pro- 
duction (NCtt"). The largest contribution to NCtt" is 
NC resonant pion production, in which the neutrino in- 
teracts with a single nucleon in the target nucleus and 
excites it to a baryon resonance; the resonant decay pro- 
duces a pion and a nucleon. NC 7r° production is the 
largest i/^-induccd background in neutrino experiments 
searching for i^^ -^ v^ oscillations. NCtt^ events cannot 
be distinguished from v^ signal events when, for example, 
one of the two photons associated with tt" — >■ 77 is not 
detected. 

Both MiniBooNE's and SciBooNE's previous measure- 
ments of NC coherent pion production were performed 
using only emitted tt'^ kinematics. However, in addition 
to the 7r° kinematics, the absence of a recoil nucleon is a 
clear and less model-dependent feature of coherent pion 
production. In SciBooNE, detection of the recoil nucleon 
is possible using the fully active and fine-grained vertex 
detector, SciBar. 

In this paper, we report a measurement of NC coherent 
tt" production using a new analysis method in which the 
lack of recoil nucleons is used to extract the fraction of co- 
herent pious within the inclusive 7r° dataset. SciBooNE's 
full neutrino data set, corresponding to 0.99 x 10^° pro- 
tons on target, is used. To simulate coherent tt produc- 
tion, the Rein and Sehgal model |6|, including lepton 
mass corrections [3], is used. The axial vector mass Ma 
and the nuclear radius parameter Rq used in the model 
are set to 1.0 GeV/c^ and 1.0 fm, respectively. These are 
the same values used in previous SciBooNE papers [3, Q ■ 
This paper updates our previous resultjOl, so, not only 
the coherent tt production model but all simulations and 
the experimental setup used in this analysis are the same 
as previously described. 



II. NC 7r° EVENT SELECTIONS 

The SciBooNE detector is comprised of three subsys- 
tems: a scintillating bar neutrino vertex detector called 
SciBar, an electromagnetic calorimeter, and a muon 



range detector. We use SciBar as the neutrino target as 
well as the particle tracker for this analysis. SciBar con- 
sists of 14336 polystyrene (CgHg) scintillator bars. The 
scintillators are arranged vertically and horizontally to 
construct a 3 x 3 x 1.7 m'^ volume with a total mass of 
15 tons. 

As shown in Fig. [H NC 7r° production is observed as 
two isolated tracks in SciBar due to two gamma rays, 
coming from the decay of the tt", converted into two e^e~ 
pairs. The background events stem from sources both 
internal and external to SciBar. Internal backgrounds 
are neutrino interactions other than NCtt*^ (mainly CC) 
within SciBar. External backgrounds come from neutrino 
interactions in the material outside of the detector vol- 
ume (dirt background events) as well as cosmic rays. To 
reduce these background events, several event selections 
are performed before extracting coherent 7r°s. 

To reject CC background events, identifying muons 
is paramount. We identify the following three types of 
tracks as muons and reject such events, (i) tracks escap- 
ing from the side of SciBar, (ii) tracks stopping in SciBar 
with delayed timing hits due to the decay electrons and 
(iii) tracks penetrating the electromagnetic calorimeter 
located at the downstream of SciBar. For the dirt back- 




FIG. 1: Event display of a typical NCtt" event candidate 
in SciBooNE data. The neutrino beam runs from left to 
right in this figure, encountering SciBar, the electromagnetic 
calorimeter (EC) and the muon range detector (MRD), in that 
order. The circles on SciBar indicate ADC hits for which the 
area of the circle is proportional to the energy deposition in 
that channel. This event display shows the electromagnetic 
shower tracks from the pair conversions of the two tt" decay 
photons. 



ground rejection, we use the upstream part of SciBar as a 
charged particle veto and require that the reconstructed 
vertices of both gamma ray candidates be in SciBar. Fi- 
nally, we require that the reconstructed invariant mass 
of two gamma ray candidates be close to the tt^ mass. 
All selections are identical to those used in the previous 
analysis Q and are described in detail there. 

After event selection, 657 events remain. Subtract- 
ing the estimated background of 240 events (202 internal 
and 38 external) yields 417 signal events. The MC ex- 
pectation is 368 events. The numbers and distributions 
obtained by the MC simulation are normalized with the 
CC data sample 0. The purity of NC 7r° production 
after all event selections is estimated to be 61%. The 
efficiency for NCtt^ production is estimated to be 5.3%. 
The efficiency for NC coherent Tr^production, incoherent 
7r° production^ with recoil neutron and with recoil proton 
are estimated to be 7.6%, 6.2% and 4.5%^. 



III. COHERENT tt" EVENT SELECTION 



incoherent tt" events is reduced from 71% in the sample 
with vertex activity to 35% in the sample without vertex 
activity. 
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In NC coherent pion production, there is no recoil nu- 
cleon in the final state since the 7r° is produced by the 
neutrino interacting with the whole nucleus. Conversely, 
a recoiling nucleon should be present in a resonant pion 
event. To separate the NC coherent 7r° events from the 
NC resonant 7r° events, recoil protons in the final state 
are used. The recoil protons are detected by their large 
energy deposition near the neutrino interaction vertex, 
so-called vertex activity. We search for the maximum 
deposited energy in a scintillator strip around the recon- 
structed vertex, an area of 40 cm x 40 cm in each view. 
The choice of 40 cm (±20 cm from the reconstructed ver- 
tex) for the area is based on the vertex resolution which is 
approximately 12 cm for each direction (x,y and z). A 7r° 
at typical SciBooNE energies travels, on average, ~20 nm 
before decaying, so the reconstructed intersection of the 
gamma tracks is a good estimate of the neutrino inter- 
action vertex. Figure [2] shows the maximum deposited 
energy distribution after all selections. Most of the co- 
herent TT*^ contribution is peaked at zero while the other 
TT*^ events have high energy activity due to recoil pro- 
tons. Events with energy deposition greater than 2 MeV 
are considered to have activity at the vertex. Note that 
incoherent pion production with a neutron recoil leaves 
no vertex activity unless the neutron kicks off protons in 
the region where we search for the energy deposit. Based 
on our MC simulation, the fraction of proton recoils in all 



^ NC incoherent tt*^ production is defined as all NCvr" events except 
for coherent Tr" production. After event selections, 89% of the 
incoherent events come from resonant pion production and the 
rest come from deep inelastic scattering. 

^ High track multiplicity around the neutrino interaction vertex 
due to the proton recoil can cause mis-reconstruction of the 
event. 



FIG. 2: Vertex activity after all event selections: the con- 
tribution from NC coherent vr", incoherent NCtt" with recoil 
neutrons, incoherent NCvr" with recoil protons, internal back- 
grounds with a tt" in the final state, internal background with- 
out a TV in the final state and "dirt" background events are 
shown separately for the MC simulation. 



IV. DATA ANALYSIS 

When a neutrino interacts with the entire nucleus, the 
following relation should be satisfied: 



> R, 



(1) 



where t and R are the four-momentum transfer to the 
target nucleus from the neutrino and the radius of target 
nucleus, respectively. This means that the cross section 
decreases rapidly when \/\t\ become smaller than R. Us- 
ing Eq. [1] we can deduce 



E-^o (1 — cos O^^o ) < — 
R 



100 MeV, 



(2) 



following Ref. [l(J|. In this equation, Et^o and Ot^o are 
the tt'^ energy and direction with respect to the neutrino 
beam, respectively. From this fact, we can determine 
the fraction of coherent tt*^ production using the recon- 
structed 7r° kinematic variable E'^°'^{1 — cosfl'^cT), where 
E'^t}^ is the reconstructed 7r° energy calculated as the sum 
of the reconstructed energies of two gamma ray candi- 
dates and O'^'i? is the reconstructed 7r° direction with re- 
spect to the neutrino beam axis. 

Wc simultaneously fit two i?'''if (1 — cos9'^°o) distribu- 
tions, with and without the vertex activity, with three 



and Ri 

^0 



scale the NC coherent 
templates independently. The 



templates made by dividing the final MC sample into NC 
coherent 7r°, NC resonant 7r° and background samples. 
Two parameters, Rcoh 
7r° and NC incoherent 
background sample is fixed to the value of the MC pre- 
diction although the systematic errors on the background 
prediction are taken into account. The expected number 
of events in the i-th bin in the E^^o{l — cos6'J^'Jf) distri- 
bution is expressed as: 



AT-p = R^„h X n: 



coh 



Rinc X N" 



■N, 



BG 



The fit minimizes the expression: 



2 ^ _„, /(jyobs.jycxp) 



(3) 



(4) 



where ]\iohs{oxp) represents the observed 
(expected) number of events in all bins 

(^1 



obs{cxp) j.^obs(cxp) 



^obs(oxp)^ and /(iV°>^^ iV^'^P) is 



the Poisson likelihood to find N°^^ events when A^^'^p 
events are expected. When the systematic errors for 
each bin and their correlation expressed with covariance 



matrix Vjk {j,k = 1,2,. 
likelihood is expressed as 



.,N{= 39))'^ are given, the 
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where A is a normalization constant. The details of the 
systematic errors and the calculation of the integral are 
described in Ref. [9|. The result of the fit is: 



Rcoh = 0.96 ±0.20, 
Ri„c = 1.24 ±0.13. 



(6) 

(7) 



The E'^°a{'^ — cosS^^cT) distribution after the fitting is 
shown in Figure [H The x^ per degree of freedom (DOF), 
before the fit is 30.8/39 = 0.79, and it is 26.6/37 = 0.72 
after the fit. Figure|4]shows three contours corresponding 
to 68%, 90% and 99% confidence level. The statistical er- 
ror and all systematic errors are included in the errors of 
Rcoh and Rinc- Without the systematic errors, we obtain 
0.98±0.18(stat.) and 1.19±0.10(stat.) for Rcoh and Rinc, 
respectively. Hence, the uncertainty of the measurement 
is dominated by the statistical uncertainty. Figures [5] and 
[S]show the distributions of the reconstructed 7r° momen- 
tum and direction with and without the vertex activity 
after fitting. 



> 

CD 



in 

CM 



<D 



LU 




Data 

NC coherent 71° 

NC other n'with n 

NC other 71° with p 

Int. BG with Jt" 

Int. BG without 11° 

Dirt 



200 400 

^r(1-cos0;!°)(MeV) 



> 

<D 



in 



Data 

NC coherent 71? 

NC other ji? with n 

NC other jtf with p 

Int. BG with jf 

Int. BG without if 

Dirt 




200 400 

^J'(1-coseJ,)(MeV) 



FIG. 3: The E'^S{1 - cos S^'b") distributions after fitting with 
(top) and without (bottom) vertex activity. 



The ratio of the NC coherent 7r° production to the total 
CC cross sections from the MC prediction based on the 
Rein and Sehgal model is 1.21 x 10"^. Hence, the cross 
section ratios are measured to be: 



cr(NCcoh7r") 
a(CC) 



Rcoh X 



cr(NCcoh7rO)A/c 



(t(CC)a/c 
= Rcoh X 1.21 X 10-^ 
= (1.16 ±0.24) X 10~^ 



(8) 



where Rcoh is 0.96±0.20. The mean neutrino energy for 
NC coherent tt" events in the sample is estimated"' to 



^ The total number of bins for the two distributions is 40 and there 
is one bin without entries. We do not include the empty bin in 
the fit. 



* In the previous paper |y|, the mean neutrino energy was 1.0 GeV 
despite using the same event sample as this paper. This is due to 
a different definition of average neutrino energy. In the previous 
paper, we used mean neutrino energy of all events passing the se- 
lection cuts in the MC simulation while, in this paper, we divide 
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FIG. 4; The contours corresponding to 68%, 90% and 99% 
confidence level for the fitted values of the scaling parameters; 
the number of degrees of freedom is 2. 



be 0.8 GeV. The fractional error of this cross section ra- 
tio is 21% while the previous result's fractional error is 
60% ((0.68 ± 0.41) X 10-2). Hence, the result has been 
improved by a factor of three with the new analysis us- 
ing vertex activity. This result is 5.8 standard deviations 
above the no coherent production assumption. The mea- 
sured cross section is also consistent with the MC pre- 
diction based on the Rein and Sehgal model [Ql- The 
result is evidence of non-zero coherent pion production 
via neutral current interactions at mean neutrino energy 
0.8 GeV. 



DISCUSSION 
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FIG. 5: Reconstructed tt momentum distributions after fit- 
ting with the vertex activity (top) and without vertex activity 
(bottom). 



A. Comparison writh the CC measurement 

The SciBooNE collaboration measured the ratio of the 
CC coherent pion to total CC production as 



g(CCcoh7r+) 
a(CC) 



(0.16±0.17(stat)iJ5:27(sys)) x lO^f?) 



at 1.1 GeV [i|- According to Eq. IHandEl the ratio of 
CC coherent pion production to NC coherent production 



the selected neutrino energy distribution by the coherent cross 
section for each neutrino energy bin before calculating the aver- 
age of the distribution. The latter method matches SciBooNE's 
CC coherent result [^ 



is measured to be 

cr(CCcoh7r+) 
cr(NCcoh7rO) 



(T(CCcoh7r+) cr(NCcoh7rO) 



cr(CC) 
n 14+0.30 



/- 



ct(CC) 



(10) 



In contrast, the Rein and Sehgal model [g as well as 
many other models predict a(CCcoh7r"'")/(T(NCcoh7r") = 
2 without the lepton mass correction |ll| . Even if we 
take it into account that the neutrino energy of the CC 
measurement (1.1 GeV) is higher than that of NC mea- 
surement (0.8 GeV), the corrected ratio ends up with 
a smaller value because the cross section increases with 
neutrino energy. So far, there is no model which can 
accommodate our measurement of the CC/NC coher- 
ent pion production ratio at these energies, although the 
measurement of the ratio at higher energies (|12|. ~7 
GeV) is consistent with 2. 
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FIG. 6: The cos (f^a distributions after fitting with the vertex 
activity (top) and without vertex activity (bottom). 



hercnt pion to NC single -k^ production is found to be: 

Rcoii X g(NCcoh7r°)MC 

Rcoh X cr (NCeoh7r")MC + Rinc X cr (NCres7r°)MC 

= (17.9 ±4.1)%, (12) 

where we assume that Rinc scales the NC single resonant 
7r° production (although Rinc actually scales all incoher- 
ent -K^ production including the multi meson production) . 
In fact, single resonant tt" production is dominates the 
incoherent tt" sample, comprising 81% after event selec- 
tions. According to Eq. [TT] and [T^l the SciBooNE mea- 
surement agrees with the MiniBooNE result within un- 
certainties. It should be noted that MiniBooNE uses a 
CH2 target and includes diffractive hydrogen scattering 
in their simulation while SciBooNE uses a CH target and 
does not include diffractive hydrogen scattering in the 
simulation. However, the effect of these differences is less 
than 10%, which is much smaller than the uncertainty of 
the SciBooNE measurement (23%). 



VI. CONCLUSION 

In conclusion, we have observed NC coherent tt*^ pro- 
duction at mean neutrino energy 0.8 GeV. The ratio of 
the NC coherent t\P production to the total CC cross sec- 
tions is measured to be 1.16 xlO~^ based on the Rein and 
Sehgal model. Our measurement confirms the previous 
MiniBooNE result. The ratio of CC coherent 7r+ to NC 
coherent tt" production is calculated to be 0.14to;28 ^^" 
ing SciBooNE's previous CC coherent pion measurement 
while many models predict 2 as this ratio. We know of 
no model that can accommodate our measurement of the 
CC/NC coherent pion production ratio. 
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B. Comparison with the MiniBooNE measurement 

The MiniBooNE collaboration measured the ratio of 
NC coherent pion to NC single 7r° production to be: 



cr(NCcoh7r°) 



cr (NCcohTrO) + cr(NCres7rO) 
(19.5 ± l.l(stat) ± 2.5(sys))%, 



(11) 



below 2.0 GeV Q, where (T(NCcoh7r*') is the cross section 
for coherent t\P production and cr(NCres7r'^) is the cross 
section for exclusive NC resonant single ■n^ production. 
The qualifier "exclusive" in the latter cross section defini- 
tion by the MiniBooNE collaboration refers to a neutrino 
interaction produced in the resonant channel and with a 
single 7r° in the final state. Using our fit result (Rcoii, 
Rinc) shown in Eq. [51 for SciBooNE, the ratio of NC co- 
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